Bounds for algorithms in differential algebra^ 



Oleg Golubitsky ^'^ Marina Kondratieva ^ Marc Moreno Maza^''^ 

Alexey Ovchinnikov ^'^ 



University of Western Ontario 
Department of Computer Science 
London, Ontario, Canada N6A 5B7 

^Moscow State University 
Department of Mechanics and Mathematics 
Lemnskie gory, Moscow, Russia, 119992 

^ University of Western Ontario 
Department of Computer Science 
London, Ontario, Canada N6A 5B1 

^North Carolina State University 
Department of Mathematics 
Raleigh, NC 27695-8205, USA 



Abstract 

We consider the Rosenfeld-Grobner algorithm for computing a regular decomposition of a radical 
differential ideal generated by a set of ordinary differential polynomials in n indeterminates. 
For a set of ordinary differential polynomials F, let M{F) be the sum of maximal orders of 
differential indeterminates occurring in F. We propose a modification of the Rosenfeld-Grobner 
algorithm, in which for every intermediate polynomial system F, the bound M{F) ^ (n — 
iy.M{Fo) holds, where Fo is the initial set of generators of the radical ideal. In particular, 
the resulting regular systems satisfy the bound. Since regular ideals can be decomposed into 
characterizable components algebraically, the bound also holds for the orders of derivatives 
occurring in a characteristic decomposition of a radical differential ideal. 

We also give an algorithm for converting a characteristic decomposition of a radical differential 
ideal from one ranking into another. This algorithm performs all differentiations in the beginning 
and then uses a purely algebraic decomposition algorithm. 
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1. Introduction 



This paper is about constructive differential algebra. We study algorithms dealing 
with algebraic differential equations. Many different problems can be addressed to this 
topic. One can, for instance, test membership to a radical differential ideal, compute 
the Kolchin dimensional polynomial. The kind of algorithms we are dealing with are 
decomposition algorithms for radical differential ideals. Generally, there are two such 
algorithms, although they have variations. 

The Ritt-Kolchin algorithm computes a prime decomposition of a radical differential 
ideal, where each prime component is represented by its characteristic set. This algorithm 
is based on important results in differential algebra (see Ritt (1950); Kolchin (1973)), such 
as the Basis Theorem, the Prime Decomposition Theorem for radical differential ideals, 
the differential version of the Hilbert Theorem of Zeroes, and the Rosenfeld Lemma. It 
also relies on the solution of the so-called factorization problem: given an autoreduced 
set, determine whether the corresponding algebraic saturated ideal is prime and, if it is 
not, find two polynomials outside of the ideal whose product belongs to the ideal. 

Due to the complexity of the factorization problem, it was desirable to avoid it , 
which was done in the Rosenfeld- Grobner algorithm proposed in ( Boulier et al. . 1995f) . 



Instead of decomposing a given radical differential ideal into prime components, this al- 
gor ithm represen t s it a s an intersection of regular differential ideals, also introduced 



in ( Boulier et al. . Il99,^ : the correctness of the algorithm, in addition to the above- 



mentioned theorems, is provided by the Lazard Lemma, which states that regular ideals 
are radical. D i ffererit proofs of this lemma can be found in ( Boulier et al. , 1997 ; Morrisonl 
19991: iHuberd . I2OOOI : iBoulier et all l2006l) . 



The Rosenfeld-Grobner algorithm is the first decomposition algorithm in differential 
algebra that has been actually implemented upto our knowledge. It forms an integral 
part of the dif f alg package in the computer algebra system Maple. Updates of this 
package are available at http : //www-sop . inria. f r/ caf e/Evelyne .Hubert/di f f alg/[ 
A more efficient implementation of this algorithm in C language can be found at the 
website http : //www . lif 1 . f r/ ~boulier/BLAD/ 

Various improvements of t he R o senfe l d-Grobner algorithm have b een proposed in 
(|Boulier et all . llQQTt iHubertl . l2000l . I2003L 120041 : iBouziane erahl . l200l[ ). They all avoid 
the factorization problem and for this reason are called factorization-free methods in 
differential algebra. However, no theoretical bound for the computational complexity of 
any of these algorithms is known. 

We make the first step towards the goal of estimating this complexity: we bound the 
orders of differential polynomials appearing in the computations. The main results of this 
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work are proven only for the ordinary case. We consider the following two bounding 
problems. The first problem is to bound the orders of all intermediate polynomials and 
the output of the Roscnfcld-Grobncr algorithm. In order to obtain such a bound in 
Proposition 10, we have modified this algorithm (see Algorithms 3 and 5) a little bit. 

It would be good to have a bound that would tell us how many times we need to 
differentiate the original system in the beginning of the algorithm, so that the rest of the 
computation can be performed by a purely algebraic decomposition algo rithm. Since f or 
algebraic decomposition algorithms complexity estimates are known (see Szantdl ( 1999t )). 
such a bound would yield a complexity estimate for the differential decomposition as 
well. In this paper, however, we do not prov ide such a b ound and, moreover, conjecture 
that it would have solved the Ritt problem We leave the discovery of such 

bound and/or the proof of this conjecture for future research. 

Nevertheless, for the second type of the algorithms we are looking at in this paper 
we obtain such a bound. Namely, we can tell how many times one needs to differentiate 
elements of a given characteristic set of a characterizable differential ideal w.r.t. one differ- 
ential ranking, in order to obtain a characteristic decomposition of this ideal w.r.t. another 
ranking. In other words, we give a bound for the conversion algorithm (Algorith m 8) for 



a characterizable ideal fr om one ranking to another (see .Boulier (1999) ; Boulier et al. 



(I2OOII) : iGolubitskvl (|2004l) for other conversion algorithms applicable to prime differen- 
tial ideals). We emphasize that the input ideal does not have to be characterizable w.r.t. 
the target ranking. We show how to obtain its new characteristic decomposition by first 
differentiating the input characteristic set and then applying only algebraic operations 
(i.e., a purely algebraic decomposition algorithm). 

The paper is organized as follows. We give an introduction into differential algebra 
in Section 2. Then we describe the original Rosenfeld-Grobner algorithm in Section 3. 
Section 4 is devoted to the bound on the orders of derivatives computed by a modified 
version of the Rosenfeld-Grobner algorithm. After that, we show how to transform a 
characteristic set of a characterizable differential ideal into a characteristic decomposition 
of this ideal w.r.t. another differential ranking. We first do this for prime differential ideals 
(Section 5) and then treat the characterizable case in Section 6. 



2. Definitions and notation 

Differential algebra studies systems of polynomial partial differential equations from 
the algebraic point of view. The approach is based on the concept of differential ring 
introduced b y Ritt. Recent t utoria l s on the constructive theory of differential ideals are 
presented in iBouhed (|200lL |2006| ): iHuberfj (|2003f) : [Si^ 12002). A differential ring is a 



commutative ring with the unity endowed with a set of derivations A = {^i, . . . ,6m}-, 
which commute pairwise. The case of A = {6} is called ordinary. If R is an ordinary 
differential ring and y £ R, we denote S'^y by y^''\ 

Construct the multiplicative monoid = |(52^(52^ • • -(5^™ | fc^ o| of derivative oper- 
ators. Let Y = {yi, ...,?/„} be a set whose elements are called differential indeterminates. 
The elements of the set QY = {9y \ 9 G Q, y G Y} are called derivatives. Derivative 
operators from Q act on derivatives as 0i(02j/i) = (^1^^2)2/4 for all 9i, 92 & & and 1 ^ i ^ n. 

The ring of differential polynomials in differential indeterminates Y over a differential 
field k is a ring of commutative polynomials with coefficients in k in the infinite set of 
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variables BY (see iKolchiiJ (|l973l ): iKondratieva et alj (|l999t ); iRitd (|l950l) ). This ring is 



denoted k{yi, . . . ,?/„} or k{y}. We consider the case of chark = only. An ideal / in 
k{y} is called differential, if for all / G / and (5 G A, 5/ G /. We denote differential 
polynomials by f,g,h,... and use letters /, J, p for ideals. 

Let F C k{yi, . . . ,y„} be a set of differential polynomials. For the differential and 
radical differential ideal generated by F in k{yi, . . . , y„}, we use notations [F] and {F}, 
respectively. 

We need the notion of reduction for algo rithmic comput ations. First, we introduce a 



ranking on the set of derivatives. A ranking ( Kolchinl . ll973h is a total order > on the set 
QY satisfying the following conditions for all 9 G Q and u,v € QY: 

(1) 0u ^ u, 

(2) u^v=^eu^ ev. 

Let M be a derivative, that is, u — 6yj for a derivative operator 

and 1 ^ j ^ n. The order of u is defined as 

ordu = ordO = ki + . . . + km- 

If / is a differential polynomial, / ^ k, then ord / denotes the maximal order of derivatives 
appearing effectively in /. 

A ranking > is called orderly iff ordw > ordu implies u > v for all derivatives u and 
V. A ranking >e/ is called an elimination ranking iff yi >ei y-j implies Oiyi >ei 02yj for 
all 61,62 e e. 

Let a ranking < be fixed. The derivative Oyj of the highest rank appearing in a 
differential polynomial / e k{?/i, . . . \ k is called the leader of /. We denote the 
leader by Id / or Uf. The indeterminate yj is called the leading variable of / and denoted 
by Iv /. Represent / as a univariate polynomial in uy: 

/ = i/u^+aiufi 



ij+aiu^ +... + ad. 



The monomial is called the rank of / and is denoted by rk/. Extend the ranking 



relation on derivatives variables to ranks: uf^ > u^^ iff either ui > U2 or ui = U2 and 
di > d2. 

The polynomial 1/ is called the initial of /. Apply any J e A to /: 



Qf 

Sf = g^Suf + Sifuj + Saiuj-^ + ... + Sad 



The leader of Sf is 5uf and the initial of Sf is called the separant of /, denoted s/. If 
9 € O \ {I}, then 0f is called a proper derivative of /. Note that the initial of any proper 
derivative of / is equal to s^. 

We say that a differential polynomial / is partially reduced w. v. t. g iff no proper deriva- 
tive of Ug appears in /. A differential polynomial / is algebraically reduced w.r.t. g iff 
degug / < *i^§ug 9- A differential polynomial / is reduced w.r.t. a differential polyno- 
mial g iff / is partially and algebraically reduced w.r.t. g. Consider any subset A C 
k{?/i ,...,?;„} \ k. We say that A is autoreduced (respectively, algebraically autoreduced) 
iff each element of A is reduced (respectively, algebraically reduced) w.r.t. all the others. 
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Every autoreduced set is finite (Kolchin, 19731 Chapter I, Section 9) (but an alge- 
braically autoreduced set in a ring of diiTerential polynomials may be infinite). For au- 
toreduced sets we use capital letters A, B, C, ... and notation A = Ai, . . . , Ap to specify 
the list of the elements of A arranged in order of increasing rank. 

We denote the sets of initials and separants of elements of A by ia and sa, respectively. 
Let Ha = ia U sa. Let S' be a finite set of differential polynomials. Denote by S°° the 
multiplicative set containing 1 and generated by S. Let / be an ideal in a commutative 
ring R. The saturated ideal I : S°° is defined as {a S i? | 3s G S°° : sa £ /}. If / is a 
differential ideal then I : S°° is also a differential ideal (see Kolchin ( 19731 )). 

Consider two polynomials / and g in k{j/i, . . . , yn}- Let / be the differential ideal gen- 
erated by g. Applying a finite number of pseudo-divisions, one can compute a differential 
partial remainder fi and a differential remainder ]i of / w.r.t. g such that there exist 
s e 5*^ and h G satisfying s/ = /i and hj = fo mod / with /i and /2 partially 
reduced and reduced w.r.t. g, respectively (see iHubertI ( 2000l) for definitions and the al- 
gorithm for computing remainders). We denote by d-rem(/, A) the differential remainder 
of a polynomial / w.r.t. an autoreduced set A. 

Let A = Ai, . . . , Ar and M = Bi, . . . , Bs he (algebraically) autoreduced sets. We say 
that A has lower rank than B if 

• there exists k ^ r,s such that rkAi — rk Bi ioi 1 ^ i < k, and rkA^ < rkSfc, 

• or if r > s and rk Ai = rk Bi for 1 ^ i ^ s. 

We say that rkA = rkB iff r = s and rkyl^ ~ rkBi for 1 ^ i ^ r. 

The following notion of a characteristic set in Kolchin's sense in characteristic zero is 
crucial in our further discussions. It was first introduced by Ritt for prime differential 
ideals, and the n extended by K olchin to arbitrary differential ideals. 
Definition 1 lKolch.iA . \l97i . page 82) An autoreduced subset of the lowest rank in a 
set X C kjl"} is called a characteristic set of X. 

We call thes e sets Kolch in characteristic sets to avoid confusion with other notions, 
e.g., in iHubertI ( 2000l . 2003 ) characteristic sets are used in Kolchin's sense and in some 
other senses. A characteristic set in Kolchin's sense exists for any set X C k|y | due to 
the fa ct that every family of autoreduced sets contains one of the least rank fsee lKolchin 
(Il973l) l. 

As it is mentioned in ( Kolchinl . Il973l Lemma 8, page 82), in the case of charfc = 0, a 
set A is a characteristic set of a proper differential ideal I iff each element of / reduces 
to zero w.r.t. A. Moreover, the leaders and the correspondent degrees of these leaders of 
any two charac teristic s ets of / coincide. 

Definition 2 iHuberi . \200ui . Definition 2.6) A differential ideal I in k{yi, . . . , j/„} is 
said to be characterizable if there exists a characteristic set A of I in Kolchin's sense 
such that / = [A] : H^. We call any such character istic set A a characterizing set of I. 
Characterizable ideals are radical ( Hubertl . 200Ct Theorem 4.4). 



3. Rosenfeld-Grobner algorithm for the ordinary case 



A system of ordinary differential equations and inequ alities A — 0, H 0, where 
A,Hc k{y}, is called regular (see iBoulier et al.l (jl995( )). if A is autoreduced, H is 
partially reduced w.r.t. A, and H D H^, where is the set of initials and separants of 
elements of A (in the partial differential case it is also required that the set A is coherent, 
but in the ordinary case this condition holds for any autoreduced set A). For a regular 
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system A, H, the differen t ial id eal [A] : H°° is also called regular. Every regular ideal is 

radical fsee lBoulier et aL and, according to the Rosenfeld Lemma, / G [A] : 

if and only if the partial remainder of / w.r.t. A b elongs to the a lgebraic ide al (A) : H°°. 

The Rosenfeld- Grobner algorithm proposed in lBoulier et alj ( 1995 . 1997 ) computes a 
regular decomposition of a given radical differential ideal {-F}, i.e., a representation 



where [A^] : are regular differential ideals. 

We begin with the following version of the Rosenfeld-Gr obner algorithm. It is very 
similar to the original algorithm presented in Boulier et al.l (|l995:), except for the fact 
that we are in the ordinary case and need not deal with coherence. We also note that 
some of the regular systems computed by the version of the algorithm presented here 
may correspond to unit ideals; th i s can be checked later on by means of Grobner basis 
computations as in Boulier et al. J 19951) or via polynomial GCD computations modulo 
regular chains as in Boulier and Lemaird ( 2000l). 

Finally, we follow the suggestion given in (|Hubertl . [2003l Improvements, page 73): it is 
recommended to reduce the multiplicative set H of initials and separants. If it turns out 
that one of them reduces to zero, then the corresponding saturated component contains 
1 and therefore need not be considered. We implement these ideas in Algorithm 1. 

Given a set F of differential polynomials, the Rosenfeld-Grobner algorithm at first 
computes a characteristic set C of F, i.e., an autoreduced subset of F of the least rank. 
It may happen that IvC C Ivi^ (for example, take F = {x + y,y} w.r.t. a ranking such 
that X > y).ln. other words, inclusion Fi C F2 does not imply that for the corresponding 
characteristic sets Ci and C2, we have Ci C C2. We need the latter property, in order to 
obtain the bound, so we are going to relax the requirement th at C is autored uced. 

A subset C of k{y} \ k is called a weak d-triangular set ( Hubertl . 2003 . Definition 
3.7), if the set of its leaders IdC is autoreduced. In the ordinary case, C is a weak d- 
triangular set if and only if the leading differential indeterminates Iv/, / G C , are al l 
distinct. A partially autoreduced weak d-triangular set is called d-triangular (jHubert . 
[2003 . Definition 3.7). For a polyno mial / and a w eak d-triangular set C, the pseudo- 
remainder d-rem(/, C) is defined via (lHubertl . [2003L Algorithm 3.13). 

We will replace the reduction of F w.r.t. an autoreduced set in the Rosenfeld-Grobner 
algorithm by that w.r.t. a weak d-triangular set. We note that the version of the algorithm 
presented in (|Hubertl . [2003I . Section 6) (Algorithms 6.8, 6.10, and 6.11) also computes 
differential pseudo-remainders w.r.t. weak d-triangular sets. Since the output regular 
systems must be partially autoreduced, at the very end, partial autoreduction of the 
weak d-triangular set C via ( Huberfj . 2003 . Algorithm 6.8) is carried out. 

Alternatively, one could perform partial autoreduction every time a weak d-triangular 
set is updated. In the following section, we show how to perform this autoreduction, as 
well as computation of differential pseudo-remainders, so that the inequality 

M{F UH)!^{n- l)!Af (Fo U Hq) 

is preserved (see formula (1) below). 
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Algorithm 1 Rosenfeld-Grobner(Fo, iJo) 
Input: finite sets of differential polynomials Fq, Hq 




{A,H)eT 



T:=0, U:={{Fo,Ho)} 
while U do 

Take and remove any {F, H) gU 

C := characteristic set of F 

^:=d-rem(F\C,C)\{0} 

H:= d-rem(iJ,C)Uifc 

if F n k = and O^H then 

if F = then T := T U {(C, H)} 

elseU :=UU{{FUC,H)} 
end if 
end if 

U:=UU{{FU {h}, H)\hG He, h^]siUH} 
end while 
return T 

4. Modified Rosenfeld-Grobner algorithm 

For a set of differential polynomials F, let mi{F) be the maximal order of the differ- 
ential indeterminate yi&Y occurring in F. If yi does not occur in F, we set mi{F) = 0. 
Let 



We propose a modification of the Rosenfeld-Grobner algorithm (see Algorithm 3 below), 
in which for every intermediate system (F, C, H) G J7, the bound 



holds, where Fq = 0, i?o 7^ is the input system of equations and inequalities corre- 
sponding to the radical differential ideal {Fg} : . 

In the formula (2) we have a multiple (n — 1)!. If the number of variables is equal to 
1 or 2 it disappears. In the case of n = 2 Ritt proved the Jacobi bound for |Fo| = 2 and 
empty Ho by the direct computation and his result docs not have any multiple either. 
Consider the intuition behind the case of n = 3 by looking at a particular example. 



n 




(1) 



M(F U C U iJ) ^ (n - l)!Af (Fo U Hq) 



(2) 
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Example 3 Let Fq = x + y+z, x' with the elimination ranking x > y > z. Then nix — 1, 
•niy = vriz — and 

M{Fo) = 1+0 + = 1. 
In order to find a characteristic set of the prime differential ideal [Fq] we reduce x' w.r.t. 
X + y + z and get y' + z'. The output consists of two polynomials: 

C ^ y' + z' ,x + y + z. 

We have: mx{C) = 0, my{C) = 1, and mz(C) = 1. Hence, 

M{C) = 0+l + l = 2>l. 

But {n - 1)! = (3 - 1)! = 2! = 2 anrf 2 2 ■ 1. 

4-1. Algebraic computation of differential remainders 

The Rosenfeld-Grobner algorithm requires to compute differential pseudo-remainders 
R = d-rem(F \ C,C). If the ranking on derivatives is not orderly, the orders of some 
(non-leading) derivatives may grow as a result of the differential pseudo-reduction, so 
that we may have mi{R) > mi{F) for some i G {l,...,n}. To ensure a bound on 
mi{R), we construct a triangular se{^ B, such that the computation of the differential 
pseudo-remainders d-rem(F, C) can be replaced by the computation of algebraic pseudo- 
remainders algrem(i^, B), and, at the same time, B satisfies a bound on the orders of 
derivatives occurring in it. 

For a set B of differential polynomials and a differential indeterminate u G IvB, let 

B„ = {/gB| \yf = v}. 

Assume that B is algebraically triangular, which implies that for any non-empty subset 
A C B, elements of A of the minimal and maximal ranks are uniquely defined and 
denoted, respectively, minA and max A. Define the following two subsets of B: 

B° = {minB„ | v e IvB} 
B* = {maxB„ | v e IvB}. 

Also, for a set {mi}^^-^ of non-negative integer numbers and an arbitrary set F of 
differential polynomials, let 

P{7ni} = {/ e F I OTdy^ / ^ rrii, i = 1, . . . , k}. 

Before we prove correctness and termination of Algorithm Differentiate<§^Autoreduce, 
let us discuss it informally. The triangular set B computed by the algorithm can be 
thought of as a result of an autoreduction of a differential prolongation of the input set 
C = {Ci, . . . , Cfc}, i.e., of the set 

C = {6^C^ 1 1 s; i fc, ^ j - d,}. 

In particular, we have rkB = rkC, unless the autoreduction process cancels one of the 
initials, in which case we can show that [C] : iJ^ = (!)• 

However, if one wants to make this autoreduction completely algebraic (in order to 
control the growth of orders), one has to be careful, because in the above set C there may 



^ A set is called triangular if the leaders of its elements are distinct. 
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Algorithm 2 Differentiate(&iAutoreduce(C, {m,}) 
Input: a weak d-triangular set C = C'l, . . . ,Ck with Id C = y['^^\ ... ,ylf''\ 

and a set of non- negative integers {mi}*^^]^, rui ^ TOi(C) 
Output: set B = | 1 < i ^ fc, < j < mj - djj satisfying 
ikBj = rkCp^ 

are reduced w.r.t. C \ {Cj} 
mj(B) < m,, i = 1, . . . , fc 

m^(l) < mi(C) + I]5'=i('^i ^ dj), i = k + l,...,n 



B C [BO] C [C] C [B] : 
HmCH^ + [C], He C {H^ + [B]) : 
or {1}, if it is detected that [C] : = (1) 

1 D:=C, B:=0 

2 while © U (5B*)|„j ^ do 

3 /:=niin(D 

^ if / e ID) then 
5 /:=algrem(/,B) 
^ D:=D\{/} 
7 else 

5 /:=algrem(/,BU((5B0\{/})) 
5 end if 

10 if rk / 7^ rk/ then return {1} end if 

11 B:=BU{/} 

12 end while 

13 return B 

appear derivatives of some IdCj of order higher than those that appear in IdC, which 

cannot be cancckd by an algebraic reduction. For example, if C = {yi, 2/2+2/1}) mi = 1, 
m2 = 2, and the ranking is elimination with yi < y2, then 

C = m, 2/1, 2/2 +2/1,2/2 +2/1:^2 +2/1 I, 

and in the last two polynomials derivatives 2/", 2/"' cannot be canceled by algebraic re- 
duction w.r.t. 2/1 and y[. 

This problem is avoided by computing the elements of B in the order of increasing 
rank. If the polynomials arc added to B in this order, one only needs to reduce each new 
polynomial /, which we are going to add to B, w.r.t. the set 5B° of first-order derivatives 
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of the elements of B° (this set BP has the same leaders as C) and the set B; this will 
guarantee that / is differentially reduced w.r.t. IB° \ Biv/. 

The inclusions! C [B°] C [C] C [B] ■.H^,HmC H^ + [C], andffc C {H^ + [m]) : 
will allow us to replace reduction w.r.t. C (in the Rosenfeld-Grobner algorithm) by that 
w.r.t. B without disturbing the saturated ideal. If m^'s arc chosen as the maximal orders 
of derivatives of y^'s, 1 ^ i ^ k, appearing in the set that is being reduced w.r.t. C, then 
we can replace the differential reduction w.r.t. C by the algebraic reduction w.r.t. B. The 
orders of derivatives of y^'s appearing in the remainder then will not exceed di for the 
leading y^'s (i.e., for 1 ^ i ^ /c). For the non- leading j/j's, the orders are bounded by the 
inequality 

k 

mi{M) ^mi{C) + '^{mj -dj), i = k + l,...,n. (3) 

We will use the following two lemmas in the proof of correctness of Algorithm DifFer- 
entiate&Autoreduce. 

Lemma 4 Let C be a weak d-triangular set in the ring of differential polynomials k{y} 
with derivations A = {^i, . . . ,6m}- Assume that a ranking on the set of derivatives QY 
is fixed. Let f e k{y} be a differential polynomial with Id/ ^ 91dC, and let f -^c 9- 
Then 

• rkg < rk/ ^ i/ £ [C] : 

• rk(7 = rk/ => 3 h € such that h ■ if — ig € [C], /i • s/ — Sg S [C]. 



Proof. Let rk / = u*^, and let A= {p G OC | Idp < u}. Then for every p € A, p and ip 
are free of u. 

Since / 9 and u ^ 01dC, there exist polynomials h £ i^, Ai, . . . ,Ak € A and 
tta, • • • ) ttfe G such that 

k 

h-f = g + Y,^i^i- (4) 

i=l 

The maximal degree of u present in (4) is equal to d. Replace every occurrence of u"^ by 

a new variable v, and consider (4) as an equality between two polynomials in v, in which 
polynomials h, Ai, . . . , A^ are free of v. We have therefore; 

h-'^ = ^ +V — A- 
dv dv ^ dv 



It remains to notice that ^ = i/ and 



dg _ ] 0, vkg<vkf 
ig, rk^ = rk/ 



dv 



hence we obtain 

• rk.g < rk/ ^ ly e (A) : C [C] : H° 



rkfii = rk/ ^ h-if-igG (A) C [C], where /i G C H^. 
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Consider now (4) as an equality between two polynomials in u, in which h, Ai, . . . , Ak 
are free of u. We have therefore: 

k 

df dg ^ dui ^ 
du du ^-^ du ' 

It remains to notice that ^ = s/ and, ifrk g = rk/, ^ = Sg, hence /i-s/ — Sg e (A) C [C], 
where h&\T C H^. □ 



Remark 5 The above lemma also holds when the set of derivations A is empty, in which 
case k{y} = k[y] is a ring of algebraic polynomials, C C k[y] is a triangular set, -^c 
is the algebraic pseudo-reduction relation w.r.t. C, and [C] = (C) is an ideal in k[y]. 
Lemma 6 (Hubert, 2 008 , Lemma 6.9) Let H and K be two sets of differential polynomi- 
als, and let I be a differential ideal. If K C {H°° + I) : , then I : ^ I : {HUK)°°. 



Proof. The proof of this statement is omitted in iHubertI (j2003l ). so, for the sake of 
completeness, we provide it here. 

Clearly, / : H°° C / : (H U K)°°. To prove the inverse inclusion, take any f ^ I : 
(H U Then, by definition, there exist h G and k e if °° such that fhk e /. 

Since K C + /) : H°°, there exist hi,h2 e H°° such that khi -h2 e I. The fact 
that fhk e / implies that fhkhi G /, whence 

fhkhi - fh{khi ~ ^,2) = fhh2 £ /, 



feI:H- 



□ 



Proposition 7 Algorithm DifferentiatecS^Autoreduce is correct and terminates. 



Proof. The proof of specifications of the algorithm is divided into three parts. We will: 

• first prove the statements about the ranks of the elements of B, 

• then about their orders, and 

• finally, the inclusions. 

All these statements hold only if the condition in line 10 is never satisfied, that is, 
throughout the algorithm rk / = rk/; for now we assume this. At the end of this proof 
we will show that, if rk / 7^ rk / for some /, then the ideal [C] : must be trivial. Let 
us proceed to the three parts of the proof: 
1. We prove that the output B has the form 

B I 1 < i fc, s$ j ^ m, - djj , 

where ABf — rkC^ . For i — 1, . . . , k, let 

[ m,(ldB) -d„ e IvB 
I —1, otherwise. 
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The statement of Part 1 follows from the following invariants of the while-loop: 



- 1 ^ ji ^ rrii - di, i = l,...,k (II) 

B= Kz^fc, 0<i<ii} (12) 

rk Bf =vkCl'\ {l^i^k, O^j ^ ji) (13) 

ivBnive = (14) 

ro c c (15) 

For all / e B,5 e ID),rk/ < rkg. (16) 



One can check immediately that the above invariants hold at the beginning of 

the first iteration of the vv'hile-loop. Assume that wc arc at the beginning of some 
iteration and the invariants hold; show that they will also hold at the end of this 
iteration. 

Let / be the polynomial computed in line 3, and let y'f'^ = Id/. We have two 
cases: 

• / e D. By 15 and the fact that the leading variables of the elements of C are 

distinct, we have rk / = rkC,;. By 14 Iv / ^ IvB, whence by definition of ji wc 
have ji = —1. Since rk / = rk/, at the end of the iteration we will have = f 
with vkBf = rkCj, and ji = 0. Thus, invariants 11-13 will hold. 

Invariants 14 and 15 also continue to hold due to the assignments in lines 6 
and 11 and the fact that sets D and B do not change elsewhere throughout the 
iteration of the while-loop. The choice of / in line 3 and the assignment in line 
6 also imply that at the end of the iteration we have rk/ = rk/ < rkgr for all 
£( G E), whence invariant 16 is preserved. 

• / G ((5B*){„.}. By 12 and 13, ji > and rk/ = rkC'f*+^\ Hence, at the end of 
the iteration ji increases by one, while / with rk/ = rk/ is added to B, thus 
preserving invariants 11-13. Note also that Iv/ G IvB, whence IvB is preserved as 
well. Hence, due to the fact that D remains unchanged throughout the iteration, 
invariants 14 and 15 are preserved. The facts that rk/ = rk/ < rkg for all 
g G D (due to the choice of / in line 3) and that Iv/ ^ IvD (due to 14) implies 
preservation of 16 at the end of the iteration. 

The above also proves the termination of the algorithm: at each iteration ex- 
actly one of the ji is incremented, whence the number of iterations does not exceed 

fc 

^(mi -rfi + 1). 

2. We first show that: 

• the elements Bj of the output B are reduced w.r.t. C \ {C,} and that 

• mj(B) ^ rrii, i = 1, . . . ,k. 

This is implied by the fact that mj(rkB) ^ rrn, i — 1, . . . ,k, which is a consequence 
of Part 1, and the following invariants: 

mt (s^) ^du {l^t^ts^k,0^j ^ J,) (17) 
B{ are differentially reduced w.r.t. B° \ {B°} (18) 
We have two cases: 
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• / e D. Let us show that 

mtif) ^ dt + jt, 2/telvB. (5) 

The fact that yt E IvB, according to 13, impHes jt ^ 0. We may assume that yt is 
present in /: otherwise nit{f) = and (5) will trivially hold due to the fact that 
and jt ^ 0. According to II, two cases are possible: 

(a) ^ jt < rrit — dt- Then there exists a polynomial g £ (5B)^^.j with 

Iv g = yt- By 13, \dg = y^'''^"''^^''. Since, due to 14, yt cannot be the leading 
variable of /, yet yt is present in /, yt"^*^^^^ < Id/. Since / is an element 
of D U (SM)'^^ y of the least rank, Id / < Idg. Combining these statements, 
we obtain 

y^'(^»<ld/<ld9 = 2/('^*+^*+^\ 

which implies (5). 

(b) jt = rrit — dt- Then, due to 15 and the condition on the input C we have 
iTT'tif) ^ rrit, which yields (5). 

Inequality (5) and invariant 17 imply that the algebraic remainder / computed 
in line 5 is differentially reduced w.r.t. B'^ and satisfies 

mt{f)^dt, ytGlvB. (6) 

Note also that due to 16, B is differentially reduced w.r.t. /. Thus, invariant 18 
also holds at the end of the iteration. 

Taking into account that for all 5 G B we have rk/ = rk/ < rk^, we obtain 

mt{f)^dt, i/tGlvD. (7) 

Together inequalities (6,7) yield invariant 17 at the end of the iteration. 

• / e ((5B*){„,^}. By 17, m,,(/) < rft + 1, for t such that yt G lvB_\ {Iv/}. This 
inequality and invariant 17 imply that the algebraic remainder / computed in 
line 8 is differentially reduced w.r.t. B° \ {Bf} and satisfies mt{f ) ^ dt- Thus, we 
obtain that invariants 17 and 18 also holds at the end of the iteration. 
Finally, we prove the bound for the orders of non-leading derivatives in the 

output: 

fe 

rns(B) < ms(C) + ^(m* - dt), s = k+l,...,n. 

i=i 

This bound holds due to the following invariant: 

m,(B)^m,(C)+ ^ jt, s = k + l,...,n. (8) 

Assume that (8) holds at the beginning of an iteration, and let s G {fc + 1, . . . , n}. 
We then have: 

(a) If / G ID, no differentiations occur during the iteration and the sum remains 
unchanged, whence (8) is preserved. 

(b) If / G ((5B*){„.}, then f = Sg for some 5 GB, whence ms{f) < ms{M) + 1. 
Similarly, 

m, (BU^B°) < m^(B) + 1. 

Thus, according to line 8, nis (./) ^ ?Ti-s(B) + 1, and so at the end of the 
iteration ms(B) is increased at most by one. At the same time, as was shown 
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in Part 1, Case 2, exactly one of the ji is incremented, whereby the sum in (8) 
increases by 1. Thus, (8) is preserved. 
3. It remains to prove: 

(a) the inchisions in the specification of the algoritlim and that 

(b) whenever the algorithm outputs {1}, the ideal [C] : is trivial. 
The inclusions are implied by the following invariants: 



B C [B°] c [C] (19) 

HmCH^ + [C] (110) 

C\Dc[B]:i?]^ (Ill) 

Hc\o C {H^ + [B]) : (112) 



Assume that the invariants hold at the beginning of some iteration; show that ei- 
ther the algorithm terminates at this iteration with the output {1} or the invariants 
will hold at the end of this iteration. We have two cases: 

• / e D. Then by 15 / G C. Since / = algrem(/,B), we have / € (B U {/}). Then, 
according to 19, / £ [C]. As was shown in Part 1, Case 1, / is added to B'^ in line 
11, thus preserving 19. 

Next, due to 14, Id / ^ IdB. Thus, Lemma 4 (see also Remark 5) applies to the 
algebraic remainder computed in line 5. We conclude from it that 

rk/>rk/ ^ ife[M]:H^. (9) 

We will use this statement later to justify the output {1}, in case the condition 
in line 10 is satisfied. For now, assume that rk/ = rk/. 
Then, from Lemma 4, we also have: 

Hj C Hf • + (B). 

Since / e C, and due to invariants 19 and 110, we thus obtain C + [C]. 
This means that 110 is also preserved. 
By definition of the algebraic remainder, 

f€{Mu{f}):H^. 

Note that line 6 results in adding / to the set C \ D, and this set is not changed 
elsewhere throughout the iteration. Thus, at the end of the iteration 111 will hold. 
Finally, as yet another consequence of Lemma 4, 

Taking into account that C \ D does not change other than in line 6, we thus 
obtain 112 at the end of the iteration. 

• / £ SM*^^_^_^ y As was shown in Part 1, Case 2, B'^ remains unchanged during 
the iteration in this case. By 19, / € [B°], whence by definition of the algebraic 
remainder applied to line 8 we have 

/ e (B U U /) C [B°] . 

Thus, 19 is preserved. 

Next, according to 12 and 13, all elements of B, and. hence, all elements of 
5B, have distinct leaders. In particular, if Id/ G ld(5B°, then / G (5B°, whence 
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Id / ^ Id ((5B° \ {/}) . In addition, since / e SB, and due to 12 and 13, we have 
Id / ^ IdB. Altogether, 

Id/^ld (BU ((5B° \ {/})) . 

Thus, Lemma 4 (see also Remark 5) applies to the algebraic remainder computed 
in line 8, yielding (9). We will use this statement later to justify line 10, assuming 
for now that rk / = rk /. 

Prom Lemma 4, we also have: 

Hf C Hf ■ H^gj^o + (B U SM^) C Hf ■ + [B]. 

Since / £ i5B, and due to invariants 19 and 110 and the fact that Hf C H^, we 
thus obtain that 110 is preserved at the end of the iteration. 

Since set C \ D remains unchanged during the iteration and set B is increased, 
invariants 111 and 112 are automatically preserved. This concludes the study of 
Case 2. 

Suppose now that rk / < rk / and apply the statement (9), which has been proved 
above in both cases. According to 19, B C [C], hence 

[B] : C [C] : C [C] : {Hn U ifc)~. 

According to 110, Hm C Hg' + [C]. Thus, Lemma 6 with H = He, K = Hm, and 
I = [C] yields [C] : (Hm U Hc)°° = [C] : H^, whence [B] : C [C] : Hg'.In 
particular, keeping (9) in mind, this implies that 

1/ e [C] : H^. (10) 

Due to 15 for Case 1, or due to 110 for Case 2, we also have that 

ifGHfCH^ + iC]. (11) 

Together (10) and (11) imply [C] : = (1). This concludes the proof of correct- 
ness. 

□ 

4-2. Final algorithm and proof of the bound 

We arc ready to present a modified version of the Roscnfeki-Grobner algorithm that 
satisfies the bound. The only place where the orders of derivatives may grow is the 
pseudoreduction w.r.t. an autoreduced set C. Of course, only the orders of non-leading 
difi'erential indeterminates may grow, while the orders of the leading ones decrease as a 
result of reduction (or stay the same if the reduction turns out to be algebraic, but then 
the orders of non- leading indeterminates do not grow either). 

By associating different weights with leading and non-leading indeterminates, we will 
achieve that the weighted sum of their orders does not increase as a result of reduction. 
These weights come from the bound in the algorithm DifFerentiateiS^Autoreduce. If the 
set of leading indeterminates changes, so do the weights. However, if we estimate in 
advance the number of times the set of leading indeterminates can change throughout 
the algorithm, we can still obtain an overall bound on the orders. 

For the original Rosenfeld-Grobner algorithm, it is not that easy to carry out such 
an estimate, because some indeterminates may disappear and reappear again among the 
leading indeterminates of the characteristic set C. For example. 
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Example 8 Let F = {y + x, x^ + z}, with the elimination ranking x > y > z. 

• We choose its characteristic set as <C :— {y + z, x}. 

• The leading variables of C are {y,x}. 

• We put F_:^d-rem{F\C,C) ^ {z}- 

• Fncw := F U C = {z, y + z, x}. 

• As radical differential ideals: 

{y + z,x,x^ + z} = [z,y + z,x] : 1°° Cl {y + z,x,x'^ + z,l} . 

• The new C — {z, x} is computed from -Fnew o.nd the leading variables have changed! 



• Finally, 

{y + z, X, x'^ + z} = [z, y, x] : 1°° ^ [z, y, x] 

and we see that the leaders y and x have come back. 
Example 9 Let F = {zy, x, x^ + z}, with the elimination ranking x > y > z. 

• We choose its characteristic set as C :— {zy, x}. 

• The leading variables of C are {y,x}. 

• We put F_:^ d-rem{F\C,C) = {z}. 

• Fncw := F U C = {z, zy, x}. 

• As radical differential ideals: 

{^zy, X, + z} = [z, zy, x] : z°° fl {zy, a:, x^ + z,z} . 

• The new C = {z, x} is computed from F^ew o,nd the leading variables have also changed! 

• But the first component is trivial: 1 S [z, zy,x] : z°° . 

The first situation can be remedied by properly relaxing the requirement that C is 
autoreduced, while the second one can detected, after which further computations in 
this branch of the Rosenfeld-Grobner algorithm are not necessary. As a result, we ob- 
tain an algorithm, in which, as long as an indeterminate appears among the leading 
indeterminates of the set C, w.r.t. which we reduce, it will stay there until the end. 

As mentioned above, we are going to replace the computation of the characteristic 
set by that of a weak d-triangular subset. It is tempting to simply compute a weak d- 
triangular subset of the least rank, since this computation is inexpensive and it would 
give us the desired property that the leading indeterminates do not disappear. However, 
the termination of the algorithm is not guaranteed then. For example, take the system 
F = {x,xy} in 'k{x,y}, and let x < y. The weak d-triangular subset of F of the least 
rank is F itself. Thus, we obtain a component {x, xy} : x°° — (1) and another component 
{x, xy,ixy}. However, i^y = x, hence we arrive at the same set F that was given in the 
input, and the algorithm runs forever. 

The reason for the above behavior is that the initials of a weak d-triangular set C, 
as opposed to an autoreduced set, need not be reduced w.r.t. C. Thus by adding these 
i nitials we do not ne cessarily dec rease the rank. The solution comes fr om the idea of 
(jBoulier et all . Il997l Section 5), (iHubertl . [2003. Algorithm 6.11), and (|Hubertl . [2004 . 



Algorithm 4.1) to construct the weak d-triangular set C gradually, so that each next 
polynomial / to be added to C is reduced w.r.t. C (thus, we can also safely add the 
initial and separant of / and guarantee that the raiik decr eases). In order to be able 
to construct the set C gradually, similarly to iHubertI (|2003h . we store it as a separate 
component of the triples {F,C,H) € U. 
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The last modification that we are going to do is the replacement of the differential 
pseudo-reduction w.r.t. C by the algebraic pseudo-reduction w.r.t. B, which is computed 
from C by Algorithm Differentiate&Autoreduce. As a result, we obtain Algorithm RG- 
Bound. 

In the proof of the bound, a key role is played by the quantity Mz{F), which is defined 
for a finite set F of differential polynomials and a proper subset Z C. Y. Assume that 
Z| = fc < n. As before, for a differential indeterminate y ^ Y, my{F) denotes the highest 
order of a derivative of y occurring in F, or zero, if y does not occur in F. Then 

Mz{F):^{n~k)J2my{F)+ ^ m,(F). 

yez yeY\z 

We also recall the notation 

Proposition 10 Algorithm 3 is correct and terminates 



Proof. We prove the following invariants of the while- loop: 
. (II) {Fo} : = HiF.cmeuiP U C} : il- n n(A,H)6T[A] : 
• For all (F, C, H) G U, 

■ (12) C is d-triangular, 

• (13) ^ is reduced w.r.t. C 

• (14) He C H, 

■ (15) Let ? = I lvC|. Then, ii I < n, 

Muc{F U C U i?) s; (n - 1) . . . (n - /) ■ M{Fo U Hq), 

otherwise 

M(F U C U F) s; (n - 1)! • M(Fo U iJo)- 
The invariants hold for the initial triple (i^Oj^i^^o)- Assuming that they hold at the 
beginning of an iteration of the while loop, we will show that the invariants also take 
place at the end of the iteration. 

Let (F, C, H) be the triple taken and removed from U . Since F 7^ 0, we can compute 
an elemen t f G F of th e least rank. Then /, as an element of F, is reduced w.r.t. C. 
Applying ( Hubertl . l2003l Proposition 6.6), we have 



{FUC}:H°°^{FUC}:iHU Hf)°° n {F U {if} U C} : n {F U {sf} U C} : 

We note that, since rki^ < rk/ and rks^ < rk/, polynomials if and Sf are, respectively, 
the elements of F U {if} and F U {sf} of the least rank (and, to repeat, their ranks 
are less than the rank of the least element of F). Moreover, since in the last two triples 
{F U {if}, C, H), {F U {sf}, C, H) only the first component has changed, invariants 12-15 
are preserved for them. For the proof of invariant II, it remains to show that 

f fB"l : F = 

{FUC} : {H\JHf)°° = { (12) 

[ {FUB"} : H°°, otherwise. 

Given that C is d-triangular, the three assignments following the computation of / ensure 
that C is a weak d-triangular set of rank strictly less than C, because the polynomial / 
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Algorithms RGBound(i^o, ^^o) 
Input: finite sets of differential polynomials Fq ^ and Hq, 

and a differential ranking 
Output: a finite set T of regular systems such that 

{Fo} : = n [A] : H°° and 

{AM)eT 

M(A UH)^{n- l)\M{Fo U Hq) for (A, H) e T. 

T:=0, U:={iFo,0,Ho)} 
while U do 

Take and remove any {F, C,H) e U 
f := an element of F of the least rank 
D:={CeC \ lvC = lv/} 
G:=FiJD\{f} 
C:=C\DU {/} 

B :=Differentiate&Autoreduce (C, {my{G UCU H) \ y e\vC}) 
if B ^ {1} then 

F := algrem(G,B) \{0} 
H := algrem(iJ,B) UiJB 
if ^ n k = and 0^ H then 

if F = then T := T U { (B°, if) } 

elseU :=UU{{F,M°,H)} 
end if 
end if 
end if 

if s/ ^ k then 

U ■.= UU{{FU{sf},C,H)} 

if 1/ ^ k then U:=UU{{FU {if}, C, H)} end if 
end if 
end while 
return T 
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is reduced w.r.t. C and we throw away (from C) all its elements with leading variables 
"in conflict" with the one of /. We note that 

GUC = (FUID) \ {/}) U (C \ B) U {/} = F U C. 

Since He C H, we also have H U Hf ~ H U H^. Therefore, 

{FUC}:{HU Hf )°° = {G U C} : (ff U ffc)°° • (13) 

Next, we use the properties of the set B ensured by Algorithm Differentiate&iAutoreduce. 
Since C + [C] , applying Lemma 6 with K = H^, we obtain 

{GUC}:{HU Hc)°° = {G U C} : {H U Hm)°° . (14) 
The inclusions B C [C] and C C [B] : imply that 

{GUC} : (HUHcU Hm)°° = {G U B} : (iJ U U Hm)°° . (15) 

Using the fact that c {H^' + [B]) : (see Algorithm 2) and applying Lemma 6 
with K = Hq, we get 

{Gum} -.{HUHcU Hm)°° = {GUM]:{H\J Hm)°°. (16) 

It follows from the definition of the algebraic pseudo-remainder (algrem) that 

{GUM}:{HUHm)'^ = {FUM}:H°°. (17) 

Indeed, {GUM} : {HU Hm)°° = {FUM} : (H U Hm)'^ . Take now any / G {FUM} : 
{H U iJB)°°- There exists h e {H U Hb)°° such that ■ f € {F UM} . li h is a remainder 
of h w.r.t. B then there exists h' e with h'h-h€ (B). Hence, 

hf e {FUB} 

and 

/ e {FUl} : H°°. 

The reverse inclusion is done in a similar way. Since B C [B°] , we obtain that {F U B} : 
= {FUBO} : 

The set B" is d-triangular, its rank is equal to that of C, set H is partially reduced 
w.r.t. B'^ and contains H^, and F is reduced w.r.t. VP. Moreover, if F = 0, we obtain the 
regular system (B'^, H), which corresponds to the radical differential ideal [B°] : H°° = 
{F U B"} : H^. Thus, we have proved (12) and also have demonstrated that invariants 
12-14 hold for the triple {F,M°,H). 

Termination of the algorithm is proved as follows. At each iteration of the while- loop, 
the triple (F, C,H) eU is replaced by at most three triples (F, B^, fl") , (F U {if}, C, H), 
and (FU{s/},C,iJ). 

Define a relation -< on the set of all triples (F, C, H) satisfying 12 -14: let (F', C, H') -< 
(F, C, H) if and only if cither rk C < rk C, or C = C and the element of the least rank 
in F' is strictly less than that in F. Then -< is a lexicographic product of two well-orders, 
which is a well-order. We have shown that in the first triple we have rkB*^ < rkC; in the 
last two triples the second component C remains the same, but the elements of the least 
rank of F U {i/} and F U {s/} are strictly less than the element of F of the least rank. 
That is, each of the three triples is less than (F, C, H) w.r.t. the well-order -<. 

Therefore, all triples computed by the algorithm can be arranged in a ternary tree, 
in which (Fq, 0, Hq) is the root, and every path starting from the root is finite. Let A 
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be the maximal length of such a path. Then the number of vertices in the tree does not 
exceed EiLo"^*- Thus, the tree is finite, whence the algorithm terminates. 

Finally, we show that invariant 15 holds for the triple (F,B°,5). We assume that 
I lvC| = I. Two cases are possible: 

(1) Iv/ e IvC. Then IvC = IvC and for any finite set of polynomials K, if I < n, we 
have 

Mi^cW = M,^c{K). (18) 

(2) Iv/ ^ IvC. Then IvC = IvC U {Iv/} and | lvC| = / + 1. If / + 1 < n, we observe 
that 











= {n-l- 


1) E rny{K)+ E rn, 


AK) = 












= {n-l- 


1) E my{K) + {n-l- 


l)-mx^f{K) + 


E my{K) 




yGlvC 






= {n-l- 


1) E my{K) + {n-l- 


2) •miv/(if)+ 






veivC 







(19) 

+ (miv/(i^)+ E rnyiK)] = 

= {n-l-l) E my{K)+ E my{K) + (n-l-2)-mi,jiK)i^ 

yelvC y^lvC 

^{n-l) E rny{K)+ E my{K) + {n - I - 2) ■ Muc{K) = 

3/elvC i/01vC 

= {n-l-l)-Mx^c{K) 

(here we have used the fact that m\vf{K) < Mivc(-^))- 
If IvC < n and | lvC| = n, we simply note that 

M{K) < Mi^c{K). (20) 

Assume for simplicity that 

where k = I or k = I + 1. Since all derivatives of yj, 1 ^ i < A;, presented in F U B U ff of 
order greater than can be found among rkB, and since the elements of F and H \ Hb, 
are algebraic pseudo-remainders of G and H w.r.t. B, we have 

rrii {F UM° U H) { (21) 
I TO,(GuBUif), k <i^n. 

Also, recall that B satisfies the inequality (see (3)) 

k 

mi(B) < mi (G U C U iJ) + ^(m^- (GuCuH) - dj), k <i^n. (22) 
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Combining (21) and (22), we obtain that 

k n 

Mi^C {F UM° [J H) = {n - k)J2di + m^{Fum°[JH) ^ 

k n 

i=l i=k+l 

k n 

^ {n- k)^d^+ ^ mi{GLlCLlH) + 

i=l i=k+l 

k 

+ {n-k) J2imj {GUCUH) - dj) = 

k n 

= (n-fc)^m, (GUCUi?) + ^ m^{GUCUH) 

i=l i=k+l 

= Mi^c (GUCUH) 

and if fc = n then 

n 

M{FuM°UH)=Ydi + = M{GuCuH) 

i=l 

because rkC = rkfiO. Thus, 

Mj^C (F UM° U H) ^ Mi^c (G U CU H) , k < n 
M (F UM^ U H) ^ M (G (J CU H) , k = n. 

Now, applying (18), (19), or (20) with K = G U CU H = F U CU H, we get 



(23) 



(24) 



Mi^c (FUM^UH) < Mivc(-F UCUH), l = k<n 

Mi^c (-P" U fiO U 5") ^{n-l- l)Mivc U ffiO U 5") ^ 

< (n-Z-l)-Mivc(i='UCUi?), l<k<n 
M{FUM°[JH) (GUCUH) = 
= M{F UCUH)^ 
^Mwc(FuCuH), l<k = n 

M(FUWPUH) (GuCuH) = 

= M{FUCUH), l = k = n. 

By taking into account the fact that invariant 15 holds for the triple {F, C, H), we thus 
obtain this invariant for the triple (F,BP,H). 

To conclude the proof of the bound for the output regular systems (M''\H), we note 
that it is already given by the invariant 15 when k = n, while in case k < n we use 
inequality (20): 

M (B° UH) ^ Mi^BO (B° UH) ^{n- 1)! • M{Fo U Hq). 
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□ 



4-3. Reduction-independent algorithm 

In Algorithm 3 wc had to be very careful in the reduction process. The idea was to 
emulate differential reductions by doing enough differentiations first and then applying 
purely algebraic reduction. We take care of the orders of derivatives in the first process 
and do not need to worry about them during the second purely algebraic step. Let us 
find out why such two-step procedure was necessary. If we reduce w.r.t. an arbitrary 
d-triangular set, the result of reduction depends on the choice of the reduction path. 
Example 11 Consider the following differential chain 

C = x{x — 1), {x — l)y, z + y + tx 

with the elimination ranking t < x < y < z and the differential polynomial 

f = z' + y'- 

We can reduce f w.r.t. C in many different ways and the remainders are very different: 
(1) 

z' + y' t'x + tx' i^^Zll^ t'{2x~l)x^2t'x^ ^t'x > 



fx /i 



(2) 



z' + y' (--^)^'+"'^ {x-l)z'-x'y (x - 1) V > 

^'+y'+''^+'^\ [x-lf{y' + fx + tx') > 0=:/2. 

We see that the remainder /i depends on the variable f that is not in both /2 and 
C. So, the reason for these so different answers is that the set C has a non-invertible 
initial. Speaking informally, if C is partially autoreduced and its initials and separants are 
invertible, then the result of reduction is more or less uniquely determined. More precisely, 
one can show that all results of reduction of a polynomial w.r.t. a d-triangular set with 
invertible initials and separants lie in a fixed Notherian ring of algebraic polynomials. In 
particular, if one of the results of reduction satisfies a certain bound on the order of its 
derivatives, then any other result of reduction will satisfy this bound as well. 

Since we are not in position of reducing w.r.t. a set with invertible initials and sepa- 
rants, we are going to state precisely and prove a slightly weaker statement. Within the 
scope of this section, let us call polynomial g a differential remainder of polynomial / 
w.r.t. C, if g is reduced w.r.t. C and there exists h € such that 

hf-ge[C]:H^. 

Proposition 12 Let C be a coheren^^ d-triangular set of differential polynomials, f a 
differential polynomial, and g a differential remainder of f w.r.t. C. Let X be the set 
of derivatives present in C and g. Let g be another differential remainder of f w.r.t. C. 



® The adjective "coherent" makes the statement vaUd in presence of partial derivatives; in the ordinary 
case, it can be ignored. 
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Assume that g is not in k[X], i.e., it admits a representation g = akvJ' + . . . +ao, where 
u ^ X and ag, . . . ,ak are free of u. Then 

5 - ao e (C) : H^. 

In particular, ao is also a differential remainder of f w.r.t. C. 

Proof. Since g and g are differential remainders of / w.r.t. g, they are both reduced 
w.r.t. C, and there exist h,hG such that 

hf-ge[C]:H^, hf-ge[C]:H^. 

Consider the differential polynomial 

/ := h{hf -g)- h{hf --g) = hg-hge [C] : . 

Since C is a coherent d-triangular set, ideal [C] : is regular. The polynomial / is 
partially reduced w.r.t. C. Therefore, by the Rosenfeld Lemma / G (C) : . 
We have 

f = {h- ak)u^ + . . . + (h- aa -h- g) 
with h ■ g contributing only to ao, because it does not depend on u. Since u does not 
appear in C, the fact that / G (C) : implies that every coefficient of / belongs to 
this ideal. In particular, h ■ ak belongs to (C) : H^, whence Oj € (C) : H^, 1 < i < fc. 
Thus, 

g-aQ = akv!^ + . . . + aiw e (C) : . 

□ 

We are going to apply the above Proposition as follows. Let C and / be as in its 
statement. Suppose we know that there exists a differential remainder g of f w.r.t. C that 
satisfies a certain bound b on the order of derivatives occurring in it. We emphasize that 
we do not need to know g, the fact of its existence is sufficient. Compute any differential 
remainder g of f w.r.t. C. Then, if g does not satisfy the bound b, it must contain a 
derivative u that does not satisfy this bound. By Proposition 12, the free coefficient of g, 
when viewed as a polynomial in u, is also a differential remainder of / w.r.t. C. Replace 
g by its free coefficient; continue such replacements until g satisfies the bound b. This 
yields an efficient procedure that computes a differential remainder satisfying the bound: 

We have proved the following 
Theorem 13 Let C be a coherent d-triangular set of differential polynomials, and let 
f he a differential polynomial. Let pi ^ mi{<C), i — 1, . . . ,n. Assume that there exists a 
differential remainder of f w.r.t. C, which contains no derivatives of differential inde- 
terminate yi of order greater than pi, i = 1, . . . ,n. Let g be any differential remainder of 
f w.r.t. C. Then Truncate(5, {pi}) is a differential remainder of f w.r.t. C, in which the 
order of every differential indeterminate yi does not exceed Pi . 

We are going to modify Algorithm 3, so that there is no necessity to perform differ- 
ential pseudo-reduction in two steps, via prolongation and purely algebraic reduction. 
In the new Algorithm 5, it is assumed that procedure d-rem computes any differential 
remainder in the above sense. The key idea is the following: whenever we find a dif- 
ferential remainder w.r.t. ED that does not satisfy the expected bound b (computed by 
Algorithm 5), by Theorem 13 we can simply truncate this remainder. In order to be able 
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Algorithm 4 Truncate {,f,{pi}) 
Input: a differential polynomial f and numbers Pi ^ 
Output: truncation of f, i.e., the sum of those terms of f that 

belong to the polynomial ring i? = k[. . . , t/j, . . . , y'f'\yi^i, . . .] 
Let / = ai + . . . + Qg, where ai are differential monomials 
g:=0 

for i := 1 to q do 

if Q, e -R then g := g + ai 
end for 
return g 

to apply Theorem 13, we are going to prove the existence of a differential remainder 
satisfying b. In fact, wc know that sets B, F, and H computed in Algorithm 3 satisfy 
b; it remains to be shown that one can obtain differential remainders w.r.t. B, given 
the elements of B, F, and H. Note that we may assume rkD = rkC (at the end of the 
for- loop), since otherwise all results of truncations arc discarded by Algorithm 5. 
To justify truncations in the for-loop of Algorithm 5, we consider 

B = DifFerentiate&Autoreduce (C, {my (G) \ y £ IvC}) . 

and show that, at the beginning of each iteration, there exist S g B and h G such 
that hB is a differential remainder of C w.r.t. P. This statement is a consequence of the 
following expanded invariant of the for-loop, which we are going to prove by induction 
on the number of iterations. Let 

B<c = {/eBl ld/<ldC}, 

B = algrem(C,B<c), E = d-rem(C,]D)), and D = Truncate {E,b). Then 

h' -C -HB ep]: H^, 
B epU{D}]: H^, 
Hb C {H^ + m) : H^, 



for some h,h' G H^. 

The inductive base holds, since at the end of the first iteration we have B = E = D = 
C and B = {C}. 

For the inductive step, we have: hi ■ C — B e (B^^) for some hi G ^m^c' 
inductive assumption [B<c] C [P] : H^. Hence, 

/ii -C-S e [B] : H^. 

Also by the inductive assumption, hi G {H-^ + W\) '■ H-^. This means that there exist 
heH^,h' € such that h ■ hi - h' e [TD]. Thus, 

h' -C-h-Bem-. H^. 
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By definition of (algebraic) pseudo-remainder, we have 

BG(B<cU{C}), Ce (£;+[©]) 

By Lemma 4, taking into account the assumption rkD = rkC, we have: 

HbCHc- H^^^ + (B<c) , He C [He + p]) : H^^ . 

By Proposition 12, E E D + (D) : H^. By modifying slightly the proof of Lemma 4, 
we will show that this implies He C Hd + (D) : H^. Indeed, using the assumption 
rkD = rkC (which holds at the end of the for- loop), we obtain rkZ) = rkC = AE; 
since all leading differential indeterminates in C are distinct, this, in particular, implies 
that 

v = ldD = ldE ^IdB. 
Now let fi, - ■ ■ ,fk be any generators of the ideal (D) : H^, so that we have 

fe 

E-D = J2c^ifi- 

By viewing the above equality as one between two polynomials in v and noting that /, do 
not involve v, we immediately obtain that ib — id S (D) : H^ and Se —Sd G (D) : H^. 

Combining the above statements, we obtain the required invariants at the end of the 
iteration: 

B e (B<c U {C}) C ([D] ■.H^U{E+ [D]) : H^) C [D U {£>}] : H^ 

and 

HbCHc- i?r<c + (B<c) C {He + [e]) : H^" + [B] : H^" C {Hd + [©]) : H^". 

The truncations applied in Algorithm 5 to compute sets F and H are justified by 
showing that differential remainders of G and H U Hf w.r.t. KD that satisfy the bound 
b exist and can be similarly obtained from the elements of sets F and H computed by 
Algorithm 3. We omit these details. 

Proposition 14 Algorithm 5 is correct and satisfies the bound. 

Proof. The proof of correctness, termination, and bound for Algorithm 5 is based on 

the same invariants of the while- loop that were used for Algorithm 3. The only new step 
we make is the Truncate algorithm whose application is justified above. □ 

5. Transformation of characteristic sets of prime differential ideals 

As above, let k{y} be a ring of ordinary differential polynomials in n indeterminates 
with the derivation S. Let C be a characteristic set of a prime differential ideal / in 
k{F} w.r.t. a ranking <. We propose an algorithm that computes a characteristic set 
of / w.r.t. any other ranking <' algebraically. More precisely, using a bound on the 
orders of derivatives occurring in the canonical characteristic set D of / w.r.t. the target 
ranking, we find a sufiicient differential prolongation of C (described below), which defines 
a prime algebraic sub-ideal / in 7 containing D. After that, it remains to compute an 
algebraic characteristic set of I w.r.t. the target ranking and extract from it a differential 
characteristic set of /. 
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Algorithm 5 RGBound-Reduction-lndependent(_Fo, i?o) 
Input: finite sets of differential polynomials Fq ^ and Hq, 

and a differential ranking 

Output: a finite set T of regular systems such that 

{Fo} : i?^ = n [A] : and 

(k,H)eT 

M(A U if) < (n - l)!M(Fo U Hq) for (A, H) G T. 
T:=0, U:={{Fo,0,Ho)} 
while U do 

Take and remove any {F, C,H) e.U 

f an element of F of the least rank 

D:={CgC \ 1vC = 1v/} 

G:=FUD\{f} 

C:=C\DU{f} 

G:=GUCUH 

b := {my (G) \ y€lvC}u\ m,, (G) + J2 {^y (G) - niy (idC)) | z ^ IvC 

B := 

for C e C increasingly do 

© := © U {Truncate (d-rem (C, D) , 6)} 
if rkD = rkC then 

F := Truncate (d-rem(G, B) \ {0}, 6) 

H := Truncate (d-rem(i? U if/, D) U Ho, h) 

if F n k = and ^ if 
then C/:=J7u{F,ID,if} 
else T := T U {(D,if)} 

end if 
end if 

if s/ k then 

U ;=[/ U {(f^U {s/},C, if)} 

if 1/ ^ k then C/ := if U {(i^ U {if}, C, H)} end if 
end if 
end while 
return T 
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5.1. A bound for characteristic sets of prime differential ideals 



First, given a characteristic set C of a prime differential ideal / w.r.t. an arbitrary 
ranking <, we would like to obtain a bound on the orders of derivatives occurring in a 
characteristic set of I w.r .t. another given r ankin g <'. For < orderly and <' arbitrary, 



such a bound is given in (|Golubitskv et al.l . 120051 ). If < is not orderly, we first obtain a 
bound for the orders o f the elements of an or derly characteristic set D of /, and then 
apply the bound from ( Golubitskv et al. . 2005f ). 

Indeed, D can be computed from C with the help of the Rosenfeld-Grobner algorithm 
applied to the system Fq ^ C, Hq = He (where the initials and separants of C in He 
are taken w.r.t. <). Since / is prime, one of the regular components {A, H) computed by 
the Rosenfeld-Grobner algorithm will coincide with /, and the characteristic set of the 
corresponding regular ideal [A\ : H°° w.r.t. <' can be extracted from the lexicographic 
Grobner basis of the algebraic ideal {A) : H°° via the algorithm given in ( Boulier et al 



I1995 . Theorem 6). A more efficient algorithm, which uses the fact that the given ideal is 
pr ime and thus avoids the computation of redundant regular components, is presented 
in lBoulier et~an (|200ll ). 

Let M be the maximal order of derivatives occurring in C. The only place where the 
Rosenfeld-Grobner algorithm differentiates polynomials is the computation of differential 
pseudo-remainders. However, for an orderly ranking, the order of a polynomial cannot 
increase as a result of pseudo-reduction. Thus, the orders of derivatives occurring in the 
characteristic set D do not exceed Al. In fact, the same applies to any other characteristic 
set of / w.r.t. the same orderly ranking: the leading derivatives of all characteristic sets 
of / w.r.t. the same ranking coincide, and the orders of non-leading derivatives occurring 
in a polynomial / cannot exceed the order of the leader of / w.r.t. an orderly ranking. 

Now we will use the following 
Lemma 15 The number of elements in a characteristic set C of a prime differential 
ideal I in the ring of ordinary differential polynomials k{?/i, . . . ,y„} does not depend on 
the ranking. 

Proof. If is a differential dimensio n of P then the number of elements of C is equal to 
n — d by (jCluzeau and Hubertl . l2003l Theorem 4.11) which does not depend on a choice 



of a differential ranking. □ 



Remark 16 The above lemma does not hold in the partial differential case. For example 
(borrowed from Boulier et al. {200ii }). a characteristic set of the prime differential ideal 



[U^ - 4U, U^j^yVy - U -I- 1, Vxx-Ux] 

in k{y} with derivations A — {d/dx, d/dy} may have 3 or 4 elements, depending on 
the ranking. 

For the above example. It takes a while to compute the characteristic set of the ideal 
w.r.t. the el i minat ion ranking u > v using the Rosenfeld-Grobner algorithm in Maple (see 
Golubitskv Consider another example that requires less computational efforts. 

Example 17 Consider the following prime differential ideal: 



P=[i 



+ y -u^+u]. 
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This set of generators forms a characteristic set of P w.r.t. the elimination ranking with 
V > u. However, if we change the ranking to u > v, then the following set containing 3 
elements will be a characteristic set of P : 

Vxxyyyj 

y ' ^xxxxyy 2y ■ V xxxxy ~t~ 2y ■ Vxxxyy ~t~ '^Vxxxx ^^xxxy H~ ^xxyy^ 
lu y • Vxxxyy ~t~ 2y • Vxxxy 2?/ ■ Vxxx '^'^xx- 



Applying Lemma 15, we obtain the following bound on the order oi I fsee lGolubitskv et al 
(l2005l) l: 

ord / := ^ ord/) ^ |C| ■ max ord C. (25) 



cec 
Den 



This bound is likely to be non-optimal. It is possible that the results of iRifrO . Il950l 



Chapter VII), together with Lemma 15, imply the following bound, which is better: let 
mi ^ TO2 ^ . . . ^ m„ be the numbers my(C), y G Y, arranged in non-increasing order, 
then 



ord/ mi. 

1=1 



For this bound, which so far is a conjecture, one needs to verify that Ritt's proof holds 
for non-elimination rankings and also adapt it for ideals specified by characteristic sets, 
rather than sets of generators. 



According to (jGolubitskv et al.U2005f ). the orders of derivatives occurring in the canon- 
ical characteristic set of / w.r.t. any ranking do not exceed the order of /. Thus, the 
number 

Ml = ICI • max ord C 

cgc 

bounds the orders of derivatives occurring in the canonical characteristic set of / w.r.t. 
any (not necessarily orderly) target ranking <'. 

We note that the bound (n — 1)! • M(C) obtained in Section 4.2 is also a bound for 
the orders of derivatives occurring in the characteristic set of / w.r.t. <' computed by 
the Rosenfeld-Grobner algorithm. In fact, invariant 15 in the proof of Proposition 10, 
together with Lemma 15, yields a better bound 

(n - IV 

{n~\C\ - 1)1 

In most cases, A/2 > Mi, but in some, especially for small values of n, it may happen 
that A/2 < Ml. This again suggests that none of the two bounds is optimal. Leaving the 
important problem of obtaining an optimal bound for future research, we summarize the 
bounds obtained so far in the following 

Lemma 18 Let <C be a characteristic set of an ordinary prime differential ideal I w.r.t. 
a ranking < . Then ord / and the orders of derivatives occurring in the canonical charac- 
teristic set of I w.r.t. another ranking <' do not exceed 

/ (n-lV 
A/c min(A/i, A//2) = min MC| • m|xordC, _ |(^| ' M{C) 
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5.2. Differential prolongation: the prime case 



Assume that ld< C = . . . |. Let nii — Mc, 1 ^ i ^ k. Compute the set 

A — DifFerentiate&Autoreduce (C, {m,j}*L J 

(for the algorithm DifFerentiatecS^Autoreduce, see Section 4.1 above). 

Let D be the canonical characteristic set of / w.r.t. <'. Every polynomial in D, as an 
element of /, reduces w.r.t. C and < to zero. Since the orders of derivatives occurring in 
D do not exceed Mc, every polynomial in D algebraically reduces to zero w.r.t. A. That 
is, D C (A) : H^. 

The algebraic ideal / = (A) : is equal to the intersection of / with the ring R = 
k [QY \ 8 ld< C U ld< A] . Indeed, A C R. Vice versa, every element oi IDR algebraically 
reduces w.r.t. A to zero and therefore belongs to (A) : H^. 

Since / is prime, so is /. Applying one of the existing efficient algorithms (for instance, 



Boulier et al. ( 2001 ) or Dahan et al. ( 20061 )) to the set A, we compute the canonical 



see 

algebraic characteristic set B of / w.r.t. the target ranking <'. We know that the algebraic 
ideal I contains the canonical characteristic set D of the differential ideal / w.r.t. <'. In 
the following section, we will show that, in fact, ID) C B. 

5.3. Extracting a differential characteristic set 

The following two lemmas hold in the partial differential case. We assume that a 
ranking is fixed. 

Lemma 19 Let k{y} be a ring of partial differential polynomials, and let K be an 
arbitrary subset of k{y} \ k. 

Let C be a differential characteristic set of K and A an algebraic characteristic set of 
K. Let T be a weak d-triangular subset of A of the least rank. Then rkT < rkC 



Proof. Suppose that a polynomial / G C is differentially reduced w.r.t. T. Then, since 
T is a weak d-triangular subset of A of the least rank, / is algebraically reduced w.r.t. 
A. Due to the fact that A is an algebraic characteristic set of K, we have / = 0, con- 
tradiction. Thus, no element of C is differentially reduced w.r.t. T, which implies that 
rkT<rkC. □ 



Lemma 20 Let L be a prime differential ideal, let C be the canonical characteristic set 
of L , and let J — /nk[y], where V C QY , be an algebraic ideal containing C. Then the 
canonical algebraic characteristic set 13 of J contains C,' more precisely, C is the weak 
d-triangular subset of D of the least rank. 



Proof. Since D is triangular, its weak d-triangular subset of the least rank is unique. 
Let T be the weak d-triangular subset of D of the least rank. 

Since ID is an algebraic characteristic set of the prime ideal J, we have Hn D J = 0. 
Moreover, Ho C k[y], therefore Ho D L — and, hence, n / = 0. Since T C / and / 
is prime, this implies 

[T] -.H^cL. (26) 

Let 

A={d-rem(/,T\{/}) | / G T}. 
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Algorithm 6 Convert_Prime (C, <, <') 
Input: a prime differential ideal P = [C] : c k{yi, . . . , y„} 

with a characteristic set C w.r.t. the input ranking < 

with leading variables yi, ■ ■ ■ ,yk cmd 

a target ranking <'. 
Output: canonical characteristic set of P w.r.t. <'. 



rui := Mc, 1 < i < A; 

A := Differentiate&Autoreduce (C, {mi}-Li) 

© := Canonical_Algebraic.CharSet ((A) : H^, <') 

return minimal d-triangular subset (D, <') 

We have A c [T] c /; we will show that set A is differentially autoreduced and rkA = 



First, show that rk A = rkT. Indeed, suppose that for some / G T and g = d-rem(/, T\ 
{/})> we have rkt; < rk/. Since T is weak d-triangular, Id/ ^ 01d(T \ {/}). Thus, 
Lemma 4 applies and tells us that i/ S [T] : . Hence, according to (26), if € I. This 
contradicts with the fact that Hj 01 — 0. 

Now, since g is reduced w.r.t. T \ {/}, rkg = rk /, and rkA = rkT, g is also reduced 
w.r.t. A \ {g}. That is, set A is autoreduced. 

By Lemma 19, rkT < rkC Therefore, rkA < rkC Since A is an autoreduced subset 
of /, while C is an autoreduced subset of / of the least rank, we have rk A > rkC Thus, 
rkA = rkT = rkC. 

Let ID) = (D \ T) U C. Set D is algebraically autoreduced, has the same rank as D, and 
satisfies the requirements of canonicity: for every / £ ID, the initial of / does not depend 
on the leaders of D, / is monic and has no factors in k[A/'(]D))], where N{Bi) = N{I}) = 
y \ IdD is the set of non- leaders of D (or D). Since the canonical characteristic set is 
unique, we have D = B) and C = T. This concludes the proof. □ 

Rciturning to the notation from the previous section and applying the above lemma, 
we obtain that the canonical characteristic set D of / is equal to the weak d-triangular 
subset of B of the least rank w.r.t. <'. This concludes the computation of the canonical 
characteristic set of / w.r.t. the target ranking, which we summarize in Algorithm 6. 

6. Transformation of characteristic decompositions of radical differential 



We generalize the algebraic method for transforming characteristic sets of a prime 
differential ideal from one ranking to another to the case of a characterizable differential 
ideal. Since an ideal characterizable w.r.t. one ranking may not be characterizable w.r.t. 
another, we need to reformulate the problem: given a characterizable differential ideal / 




rkT. 



ideals 
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with a characteristic set C w.r.t. a ranking <, compute a characteristic decomposition of 
/ w.r.t. another ranking <' algebraically. By analogy with the prime case, an algebraic 
computation here means finding a sufficient differential prolongation of C, which defines 
a characterizable algebraic sub-ideal / in /, such that a differential characteristic decom- 
position of / w.r.t. <' can be extracted from an algebraic characteristic decomposition 
of /w.r.t. <'. 

We note that, given a characteristic decomposition of a radical differential ideal w.r.t. 
one ranking, we can obtain its characteristic decomposition w.r.t. another ranking alge- 
braically by solving the above problem for each characterizable component. 

All results of this section hold in the partial differential case, except for the boimd in 
Section 6.2, which so far is known only for the ordinary case. 

6.1. Differential prolongation 

Definition 21 Let F be a subset in a ring kjK} of partial differential polynomials with 
a set of derivations A. A set G C QF is called a differential prolongation of F , if 
F <Z G and the complement of G, QF\G, is invariant w.r.t. differentiation, i.e., for all 
f eeF\G andd e A, 6f eeF\G. 

A particular case of a differential prolongation of a weak d-triangul ar set F is 
itself. If F = C is autoreduced an d coherent then, according to (Kolchinl. I1973L Lemma 



6, page 137) and ( Hubert . 2000l Lemma 6.1 and Theorem 6.2), the differential ideal 



/ = [C] : is prime, respectively characterizable iff the algebraic ideal J = (C) : 
is prime, respectively characterizable. The ideal J can be considered either as an algebraic 
ideal in the ring of differential polynomials k{y} or as an ideal in the polynomial subring 
k[Zc], where Zc = LU N, L = \dC, N = QY \ QL, since the fact that C is autoreduced 
implies C C k[Zc]. The Rosenfeld Lemma states that 

[C]:i/^nk[Zc] = (C):ff^, 

where the latter ideal is considered in k[Zc]. Moreover, a set D is a differential charac- 
teristic set of / iff B is an algebraic characteristic set of J (if the latter is considered 
in k[Zc], otherwise we need to impose an additional requirement that D is differentially 
autoreduced). In particular, the canonical characteristic sets of / and J (differential and 
algebraic, respectively) coincide (for this statement, it does not matter in which ring to 
consider J, since the canonical characteristic set of an ideal is the same regardless of the 
ring in which the ideal is considered). 

Now, if we consider a differential prolongation D of C and the corresponding polyno- 
mial subring k[Zn], where Zn = LU N, L = km, N ^ BY \eL ^ BY \eL, then B is 
not necessarily a subset of k[ZD]: 

Example 22 Let C = y' , x + y with the elimination ranking y < x and a prolongation 

M) = y,x + y, x+y,x +y . 

Then 

L = y', X, x', x", N = y. 
Hence, we have that x" -\- y" ^ k[ZD]. Also, 

[C]:H^nk[Zji]^{y', x + y, x' , x") 

and x" ^ (©) : . 
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Therefore, we need to distinguish between two ideals lo ■= (©) : in k{Y} and 
Id :=/nk[ZD] ink[ZD]. 

The algebraic ideal In depends only on the set of leaders L of the differential prolonga- 
tion of C. In other words, for any characterizing set C of / and its differential prolongation 
D with IdD = IdD = L, we have /p = Id- We call := Id a prolongation ideal of the 
ideal I. 

Next, we study the properties of the prolongation ideals. The following lemma gives 
a criterion for a prolongation ideal to be prime or characterizable. 
Lemma 23 Let C be a coherent autoreduced set, and let H be a differential prolongation 

of C Then the differential ideal 1 ^ J = [C] : is prime, respectively characterizable, 
iff the corresponding prolongation ideal Id is prime, respectively characterizable. 

Proof. If / is prime then its restriction / n k[ZD] = lo is also prime. If lo is prime than 
its restriction Id fl k[Zc] = (C) : is prime and, thus, / is prime. 

Let 7 be a characterizable differential ideal. We will show that set A given by for- 
mula (27) characterizes the prolongation ideal lo- We have lo C (A) : H^. Indeed, by 
Lemma 4, sets A and D have the same ranks, whence they have the same sets of reduced 
polynomials. In particular, since D is a differential prolongation of the characteristic set 
C, the ideal In has no non-zero polynomials reduced w.r.t. D, and hence w.r.t. A. 

Now note that (A) : C / and A C k[Zo]. Hence, Id = (A) : and A is a 
characteristic set of lo. Thus, lo is characterizable. 

Since C C \i[Zo] and {C) : = I n k[Zc], we have 

{£) : = {A) : nk[Zc]. 

□ 

The next lemma establishes a relation between the characteristic sets of a character- 
izable differential ideal / and the algebraic characteristic sets of its prolongation ideals. 

Lemma 24 Let C be a characteristic set of the differential ideal 1 ^ / = [C] : , let 
L be a differential prolongation of L = \d C, and let Ii be the corresponding prolongation 
ideal. 

Then a characterizing set A of Ii can be obtained from C as 

A := {algrem(/, B \ {/}) | / G B, Id / e L}, (27) 

where B is any triangular subset of Q<C satisfying IdB = IdGC. 

Vice versa, given a characterizing set A of Ii, let T be a weak d-triangular subset of 
A of the least rank. If T is differentially autoreduced, then it is a characterizing set of 
I. In particular, if A is the canonical characteristic set of Ii, then T is the canonical 
characteristic set of I. 

Proof. Since / is characterizable, Ii is also characterizable by Lemma 23 and A is its 
characteristic set. The other way follows from Lemma 19. □ 

In the ordinary case, the triangular set B considered in the above lemma is unique. 
Moreover, set A can be equivalently obtained as 

A := Differentiate&Autoreduce(C, {toi}). 
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where the numbers {m^} are the maximal orders of derivatives of the leading differential 
indeterminates of C occm'ring in the prolongation L. It is preferable to compute A in 
this way, because DifFerentiatecS^Autoreduce provides a bound on the orders of non-leading 
derivatives occurring in A, which can be used for establishing complexity estimates for 
the entire transformation algorithm. 

A generalization of Algorithm Differentiated Autoreduce to the partial case is an inter- 
esting open problem. Moreover, in the partial case, there may be uncountably infinitely 
many triangular subsets of 0C whose leaders coincide with Id 0C. Thus, not every such 
set can be enumerated by an algorithmic procedure. However, it is easy to write a proce- 
dure that would enumerate a particular subset of OC, given C; this procedure makes the 
computation of the set of algebraic pseudo-remainders algorithmic as well. If one would 
like to choose the subset B in a systematic way, we suggest to use t he ideas from the 
theory of monom ial involutive d ivisions (see lGerdt and Blinkov ( 19981 )). 



According to ([Hubert, 2003 . Theorem 4.13), there is a one-to-one correspondence 
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between the essential prime components of a characterizable differential ideal [C] : H^ 
and the minimal prime components of the corresponding algebraic ideal (C) : . The 
following lemma generalizes this result to prolongation ideals. 

Lemma 25 Let C be a characteristic set of the differential ideal / = [C] : , let L 
he a differential prolongation of L = Id C, and let Ii be the corresponding prolongation 
ideal. 

Let / = Pi n . . . n Pfc be the essential prime decomposition of L, and let {Pi)i be the 
prolongation ideals corresponding to Pi, i — 1, . . . , k. Then 

h = {Pi)i n . . . n {Pk)L 

is the minimal prime decomposition of L^ . 
Proof. Since /£ = / n k[Z£], 

/£ = (Fi n k[Z£]) n . . . n (Pfc n k[Z£]) = n . . . n {Pk)i 

is a prime decomposition of the ideal J^. Suppose that it is not minimal. Since (C) : 

i?^ = /£nk[Zc], 

(C) : = ((Pi)£ n k[Zc]) n . . . n ((Pfc)^ n k[Zc]) 

is a prime decomposition of the ideal (C) : , which is also not minimal. But the latter 
contradicts the fact that {Pi)i fl k[Zc] = PiCi k[Zc], 1 ^ i ^ fc, and 

(C) : = (Pi n k[Zc]) n . . . n (Pfc n k[Zc]) 

is the minimal prime decomposition. □ 



6.2. A bound for characteristic sets of prime components 

Let / = [C] : be a characterizable differential ideal with a characteristic set C 
w.r.t. a ranking <. Let L = ld< C, and let L be a differential prolongation of L. From the 
previous section we know that the prolongation ideal L^ is characterizable (Lemma 23) 
and its minimal prime components correspond to the essential prime components of / 
(Lemma 25). We would like to find a sufficient differential prolongation L such that the 
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minimal prime components of contain differential characteristic sets of the correspond- 
ing essential prime compone nts of I w.r.t. a ny other ranking <'. 



First of all, according to (jHubertl . l2003l Theorem 4.13), a differential characteristic 



set of an essential prime component of / coincides with an algebraic characteristic set of 
the corresponding minimal prime component of the ideal (C) : . This implies that 
every essential prime component P oi I has a characteristic set C p satisfying the bound 
niyiCp) ^ my(C) on the orders of derivatives of any differential indeterminate y Cz Y 
occurring in Cp. 

For the ordinary case, as was shown in Section 5.1, we thus have a bound Mc on the 
orders of derivatives occurring in the canonical characteristic sets of the essential prime 
components of / w.r.t. any other ranking <'. For the partial differential case, such a 
bound is not known, but let us assume that we can compute such a bound Mc also for 
the partial caseQ We need to assume that Mc ^ niy{<C) for all y gY. 

Let 

L = {eu\ueL, OTd0u!^Mc} (28) 
be the differential prolongation of L up to the order Mc- According to Lemma 25, the 
minimal prime components of contain all polynomials of the corresponding essential 
prime components of / of order less than or equal to Mc- Thus, they also contain the 
canonical characteristic sets of the corresponding essential prime components of / w.r.t. 
any other ranking <'. In what follows, we will denote the above differential prolongation 
II simply by /. Applying Lemma 24, we compute a characteristic set of / w.r.t. <. 

6.3. Algebraic hi- characteristic decomposition 

So, we have the differential ideal / which is characterizable w.r.t. the ranking < and 
would like to give a characteristic decomposition of / w.r.t. <' . We have constructed the 
prolongation algebraic ideal / which is characterizable w.r.t. < with a characteristic set 
A given by formula (27). Let 

/ = Ji n . . . n Jfc (29) 

be a bi-characteristic decomposition of / w.r.t. < and <'. That is, each component Ji, 
1 ^ I ^ /c, is an algebraic ideal characterizable w.r.t. both rankings with the canonical 
characteristic sets and w.r.t. < and <', respectively. 

Let us discuss how one can construct such a decomposition. Algorithm 7 does the 
following. Given a characterizable algebraic ideal / with the characterizing set C w.r.t. 
<s, it first computes its (possibly redundant) algebraic characteristic decomposition w.r.t. 
<t via the procedure 

Algebraic-characteristic-decomposition{<C, <s, <t). 

This procedure can be performed, for example, by applying the Triade algo rithm dMoreno MazaL 



119991) . which is implemented in the RegularChains library in Maple fsee lLemaire et al 
(j2005r )V A parallel implementati on of t his alg orithm , on a shared memory machine in 
Aldor is also in progress (see Moreno Maza and Xiel (j2006h ). 



If one of the characterizable components turns out to be equal to / (note that equality 
of characterizable algebraic ideals can be checked, e.g., by computing their Grobner 



^ Of course, Mc can be obtained by computing characteristic sets of the prime components w.r.t. the 
target ranking, but this would clearly defeat our purpose: we need a bound that can be computed from 
C relatively easily. 
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Algorithm 7 Algebraic-Bicharacteristic-Decomposition (C, <, <') 
Input: characterizing set C of a characterizable algebraic ideal I 
w.r.t. an ordering < on variables 
and another ordering <' 
Output: a finite set T = {(Ci,ID)i) | i e 3}, where 

for every i € 3, Ci and are algebraic characterizing sets 
of the same ideal li w.r.t. < and <', respectively, and 
I = Hiejli 

<s:=<, <*:=<' 
£ := {C}, T := 
while £ / do 

U :=€,€:= 

for C e J7 do 

J := (C) : w.r.t. <s 

S) :=Algebraic-characteristic-decomposition{C, <s, <t) 
if 3 D e S such that J = (D) : w.r.t. <t then 

if <,s=< then T := T U {(C, D)} else T := T U {(D, C)} 
else £ := £ U 2) 
end if 
end for 

if <,,=< then <«:=<', <t:=< else <«:=<, <t:=<' 
end while 
return T 

bases), then / is bi-characterizablc; in this case the algorithm terminates and outputs T 
consisting of a single pair (C,©) of characterizing sets of / w.r.t. < and <', respectively. 
If all characterizable components of / contain it strictly, then, for each characterizable 
component, we compute its characteristic decomposition w.r.t. < and repeat the above 
strategy. 

Correctness of the algorithm follows from the fact that, at each iteration of the while- 
loop, CUT provides a characteristic decomposition of / w.r.t. and T satisfies the 
requirements of the output. Termination follows from the Notherian property of the 
polynomial ring, i.e., that every sequence of strictly nested polynomial ideals is finite. 

We note that components J^, for which ld< A,; =^ ld< A, arc redundant, i.e., they can 
be excluded from the right-hand side of (29) without aflFecting the intersection. Indeed, 
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if / = Pi n . . . n P/ is the minimal prime decomposition of /, and Ji — Qi.i H . . . fl Qi^i^ are 
the minimal prime decompositions of J^, 1 ^ i ^ fc, then a component Ji is redundant, 
if none of Pj, 1 ^ j ^ I, can be found among Qi^t, ^ ^ t ^ li. But this is the case if 
ld< Aj 7^ ld< A, since by ( Hubert . 20031 Theorem 4.13) the characteristic sets of Pj have 
leaders ld< A, while the characteristic sets of Qj,* have leaders ld< A^. Therefore, we can 
assume that for all 1 ^ i ^ fc, ld< A^ = ld< A. 

We prove then that every minimal prime component of Ji is a minimal prime com- 
ponent of /. Indeed, every Qi^t is a prime ideal containing /. Suppose that Qi^t is not 
minimal, i.e., there is a minimal prime component Pj of / such that Pj C Qi,t- But the 
latter strict inclusion is impossible according to the following Lemma 26 and Remark 27. 
Lemma 26 Let P and Q be two prime differential ideals whose characteristic sets w.r.t. 
< have the same sets of leaders Then P <^ Q implies P = Q. 



Proof. Let Ci and C2 be these characteristic sets. We have P — [Ci] : and Q = 
[C2] : Consider the restricted ideals p = (Ci) : i7g and q = (C2) : in 

the Notherian ring k[L, A''(Ci, C2)], where A^(Ci, C 2) is the set of non-leading variables 
appearing in both Ci and C2. From ( Hubert! . 200d Theorem 3.2) it follows that both p 
and q are of dimension \N{Ci,C2)\- 

Take any / G p. It is partially reduced w.r.t. both Ci and C2 (which are coherent 
and autoreduced) and belongs to P C Q. By the Rosenfeld lemma / G q. Hence, p C q 
and they are prime and must be equal then, because their KruU dimensions are equal 
to the same number |7V(Ci ,C2)|. Hence, we have C i C Q and C2 C p C P at the same 
time. Thus, according to ( Golubitsky et al. . 20051 Theorem 9) we finally obtain that 
P = Q. □ 



Remark 27 In the above lemma, one can assume that the set of derivations is empty, 
hence the statement also holds for algebraic ideals. 

To summarize, for every bi-characterizable component Ji, there exists a subset Ti C 
{I, . . . ,1} such that 

is the minimal prime decomposition of Ji. Moreover, equality (29) implies that 

|Jt, = {i,...,o. 

6.4. Constructing differential characterizable components from the algebraic ones 

Fix any of the above algebraic bi-characterizable components J — Ji, where 1 ^ i ^ k; 
we have a set of indices T — Ti C {1,...,^} such that 

As above, let A = A^ and B = be the canonical characteristic sets of J w.r.t. < and 
<', respectively. 
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According to Lemma 25, each minimal prime component Pj of / is a prolongation 
ideal of the corresponding essential prime component Pj of /, i.e., 



Pj = Pj n k[L U N] 



where / — Hj^i the essential prime decomposition of /. Since IB is a characterizing 
set of J w.r.t. <', the initials and separants of B w.r.t. <' are not zero-divisors modulo 
J, i.e., they do not belong to the minimal prime components Pj, j £ T. Since B, as well 
as iJfl, is a subset of k[Z U A^], we have therefore Hr Ci Pj — 0, j £T. 

Let T C B be the weak d-triangular subset of B of the least rank w.r.t. <'. Since 
Hj C Hm, we also have Hj n Pj = 0, j e T. Thus, [T] : C Pj, j e T. In particular, 
this implies that [T] : HJ^ ^ (1). 

Let D be the result of differential autoreduction of T w.r.t. <', i.e., 

D = {d-rem(/,T\{/}) | / £ T}. 

Set D is differentially autoreduced. By definition of differential remainder, D C [T]. By 
Lemma 4, since [T] : ^ (1), we have rk</ ID) = rk</ T and, moreover, Hjj, C + [T]. 
Therefore, 

m-.H^ C[T]:H^ cPj, jeT. (30) 
We will show that D is a characteristic set of the ideal [D] : w.r.t. <' by proving 
that every polynomial in the intersection PljeT reduces w.r.t. D to zero. Given (30), 
this will also imply that 

[B] ■.H^=f]Pj. (31) 

Take any polynomial / e rijeT^i' ^'^'^ / ~ d-rem(/, D), where the pseudo- 
remainder is computed w.r.t. <'. Since D C [T] C Pj, j G T, we have / G Hjer Pj- 

Let ¥j be the canonical characteristic set of Pj w.r.t. <', and let Fj be the canonical 
algebraic characteristic set of the corresponding prolongation ideal Pj . We have shown in 
Section 6.2 that Pj contains ¥j. Thus, from Lemma 20 it follows that ¥j is the weak d- 
triangular subset of F^- of the least ra nk w.r.t . <'. O n the other hand, since Pj is a minimal 
prime component of J, according to (Hubert, 20031 Theorem 4.13), ld<' ¥j = ld<' B. This 



implies that ld<' ¥j = ld<' T = ld</ D. That is, the fact that / is reduced w.r.t. D implies 
that it is partially reduced w.r.t. F^. 
By the Rosenfeld Lemma, 

/G(F,):i?r7c(F,):i7- = P„ jeT 

i.e., f G J. Now, the fact that / is reduced w.r.t. D implies that it is algebraically reduced 
w.r.t. B. Since the latter is a characteristic set of J, we obtain f — and the required 
equality (31). 

Now we see that the ideal [D] : is characterizable w.r.t. <'. The canonical char- 
acteristic set of this ideal w.r.t. <' is contained in each minimal prime component of the 
ideal (ID) : H^, therefore it is also contained in every Pj, j G T, and hence in J. The 
ideal J is contained in [D] : H-^. Thus, by Lemma 20, the canonical characteristic set of 
[ID] : is equal to the weak d-triangular subset of B of the least rank w.r.t. <'. That 
is, we have 

D = T 
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Algorithm 8 Convert_Characterizable (C, <, <') 
Input: set C which characterizes the ideal [C] : w.r.t. the input ranking < 

and has leading variables yi, ■ ■ ■ ,yk and a target ranking <' . 
Output: characteristic decomposition of [C\ : w.r.t. <'. 

Mc := min (|C| • m|xordC, („L"|c|2i), • M(C)) 

m, := Mc, 1 < i < fc 

A := Differentiate&Autoreduce (C, {mijf^i) 

D := Bi-characterizable_CanonicaLDecomposition ((A) : H^, <, <') 
€ := {minimal d - triangular subset (©, <')\B& D} 
return £ 

which is (w.r.t. the ranking <') the canonical characteristic set of the characterizable 
differential ideal 

[D] : H^. 

6. 5. The final characteristic decomposition 

In the previous section, we have shown that for each bi-characterizable component 
Ji, 1 < i < I, of I with the canonical characteristic set ffij w.r.t. <', if Dj is the weak 
d-triangular subset of of the least rank, then it is the canonical characteristic set of 
the ideal [Di] : H^.. We have also shown that 

m ■.Hs:= f]Pj. 

jeTi 

Thus, since lJi=i Ti = {1, . . . , I}, the following intersection 

i=l 

is a characteristic decomposition of / = Pifl. . .(iPi w.r.t. <'. This concludes the algebraic 
computation of a characteristic decomposition of / w.r.t. the target ranking, which we 
summarize in the Algorithm 8. 

Now, in order to convert a characteristic decomposition 

/=f][Q]:H^ 

i=l 

of a radical differential ideal / w.r.t. < to a ranking <', one just applies Algorithm 8 to 
each characterizable component [Cj] : and then collects all the results together in a 
single intersection. 

7. Conclusions 

By estimating the orders of derivatives, we have shown that, given a set of ordi- 
nary differential polynomials specifying a radical differential ideal /, one can construct 
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a Nothcrian ring of algebraic polynomials, in which the computation of a characteristic 
decomposition of / is actually performed. This does not mean that the computation is 
completely algebraic: differentiations are allowed, but they never lead out of the con- 
structed algebraic ring. 

For the problem of converting a characteristic decomposition of a radical differential 
ideal from one ranking to another, we have proposed an algorithm, which first differen- 
tiates the input polynomials sufficiently many times, and then performs the conversion 
completely algebraically, without using differentiation at all. The algorithm is applicable 
in the partial differential case, but the bound for the number of differentiations of the 
input polynomials is given for the ordinary case only. 

We conjecture that, if one can solve the first problem of computing a characteristic 
decomposition of a radical differential ideal from generators completely algebraically, 
i.e., by an algorithm that first differentiates the input polynomials sufficiently many 
times, and then computes the decomposition without using differentiations, then one 
can also solve the Ritt problem of computing an irredundant prime (or characteristic) 
decomposition of a radical differential ideal. 
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